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ABSTRACT 

The binary YY Gem shows many interesting properties, one of which is the periodicity 
in its flaring rate. The period, which is about 48 ±3 min, was ever interpreted in terms 
of the oscillation of a filament. In this paper, we propose a new model to explain this 
phenomenon by means of 2.5-dimensional MHD numerical simulations. It is found that 
magnetic reconnection is induced as the coronal loops rooted on both stars inflate and 
approach each other, which is driven by the differential stellar rotation. The magnetic 
reconnection is modulated by fast-mode magnetoacoustic waves which are trapped 
between the surfaces of the two stars, so that the reconnection rate presents a periodic 
behaviour. With the typical parameters for the binary system, the observed period 
can be reproduced. We also derive an empirical formula to relate the period of the 
flaring rate to the coronal temperature and density, as well as the magnetic field. 

Key words: binaries: close - stars: flare - stars: oscillations - waves - methods: 
numerical 



1 INTRODUCTION 

The binary syste m YY Gem has been s tudied for many years 
since discovered (|Adams fc Jovlll920r ). It is a double- lined 
spectroscopic eclipsing binary system which contains two 
dMle late-type stars with masses and radii a lmost identical, 
and t he orbital perio d being 0.8142822 d (|jov fe Sanfordl 
1 19261 : Ivan Genu fl93ll ) . It is of great importance to study 
the basic relations of such rare lat e-type binaries l|Kronl 
1 19521 : iButler. Doyle fc Buddindll996h. Based on some con- 



siderations in iDworakl l|l975l h1 Brancewicz fc Dworakl (l980) 
gave the parameters of YY Gem: Mi = M 2 = 0.57 M Q , 
Ri = 7?2 =0.6 R,0, the separation of the two components 

a = 3.83 R Q . 

IChabrier fc Baraffd (1995) demonstrated that the depth 
of each stellar radiative core is ~ 70 per cent of its 
radius, and the thickness of the convective zone is ~ 
30 per cent. Strong convective motions under the sur- 
faces cause large-scale star spots and huge magnetic field 
structures. Flare activity on YY Gem was reported for 
the first time by Moff ett fc Boppl l|l97ll ). Subsequent in- 
vestigations have sh own that it is one of the most ac- 
tive flaring binar ies (|MoffettJ 1 19741 : iDovle fc Butlerl 1 19851 : 
iDovle et al.lll990l ). Their flaring activities have been stud- 
ied in multi-wavelengths . Usin g the Very Large Array, 
Ijacksqn. Kundu fc White] (1 19891 ) presented radio observa- 
tions: Istelzer et al.l (1200*21 ) studied the simultaneous X-ray 
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spectroscopy of YY Gem with Chandra and XMM-Newton 
satellites; While, far UV observation s of YY Gem were re- 
ported by ISaar fc Bookbinderl l|2003l) . The flarings on YY 
Gem exhibit UV and X -ray emissions whi ch are stronger 
than those on the Sun (|Haisch et al.l fl990). For example, 
iTsikoudi fc Kellettl (|2000l ) observed two large flares, and the 
integrated X-ray luminosity was estimated to be about 6-8 
x 10 33 erg. They also calculated the ratio of X-ray and bolo- 
metric luminosities Lx/LboI, which indicated strong mag- 
netic activities and 'hot' coronal components. 

For close binaries with similar masses, the density scale- 
height is quite large, i.e., the density does not decrease 
with height from the stellar surfaces as rapidly as in the 
solar corona. Therefore, the hot plasma between the two 
stars could emit strong X-ray light. For example, observa- 
tions revealed large-scale strong X-ray so urces between the 
two stars of the RS CV n binary AR Lac (|Siarkowskilll992l : 
ISiarkowski et al.l I 19961) as well as of the RS CVn binary 
TY P yx <|Culhane et al.lll990l : [Pres. Siarkowski fc Svlwesterl 
Il995l ). At the same time, it is probable th at the two stars are 
magne t ically connected, as illustr ated bv lUchida fc Sakurail 
(1983). lUchida fc Sakurail (|l985l ) proposed that the inter- 
stellar activities are intimately affected by the differential 
rotations of both stars. So, there may exist some interesting 
phenomena in such a system. On e of them is the p eriodic 
flaring rate, which was reported bv lDovle et"aT1 (|l99ol ). Their 
observation in the [/-band showed four flares separated in 
succession with a periodicity of 48 ± 3 min during a total 
observing time of 408 min. The duration of each flare varies 
from ~ 20 min to ~ 40 min, and the time-averaged flare 
luminosity in the [/-band is 1.25 x 10 28 erg s~ . They in- 
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terpreted this periodicity in terms of filament oscillations. 
In this paper, we attempt to propose an alternative expla- 
nation for the periodicity. 

Flares on the Sun and dKe-dMe stars are both thought 
to result from magnetic reconnection. On the Sun, periodic 
behaviours have also been reported in various reconnection- 
associated phenomena, e.g., the quasi-periodic modula- 
tion of flaring emission l|Nakariakov et all l2006h and the 
repetitive appearance of transition region explosive events 
l|Chen fc Pr iest 2006). In the former case, magnetic recon- 
nection is modulated by magnetoacoustic waves from a 
nearby oscillating coronal loop, whereas in the latter, mag- 
netic reconnection is modulated by p-mode waves. 

To our knowledge, no efforts in MHD numerical simula- 
tions have been devoted to the investigation of flaring phe- 
nomenon in binary systems, although hydrodynamic simula- 
tions have been done rega rding flares in M dwarf stars (e.g. 
ICheng fc Pallavicinilll99ll). Using potential field extrapola- 



Uchida fc Sakurail (|l985h demonstrated that the evolv- 



tion. 

ing surface magnetic field can result in magnetic reconnec- 
tion in the interstellar space. Here, we perform MHD numer- 
ical simulations in order to investigate the magnetic recon- 
nection in a close binary system, with the purpose to clarify 
which process is likely to be r esponsible for the p eriodicity of 
the flaring rate as reported bv lDovle et al.l l|l990l ). The paper 
is organised as follows. Section 2 describes the basic model, 
MHD equations, initial and boundary conditions. Section 
3 gives the numerically calculated results. Discussions are 
presented in section 4, which is followed by conclusions in 
section 5. 



2 NUMERICAL METHOD 
2.1 Problem setup 

The basic model used in our work is shown in Fig. 1. For 
ordinary binaries, the member stars are usually of different 
types. As usual, we call the one with larger mass the pri- 
mary, marked as Mi, the other the secondary, marked as 
M2 (Mi > M2). Owing to the difference in mass, the point 
of zero gravity and the centre of mass are located at different 
positions. The former is closer to the secondary (M2) and the 
latter is closer to the primary (Mi). In the case of YY Gem, 
we ta ke the parameters given by iBrancewicz fc Dworakl 
(1980), i.e., Mi = A/2, so that the two points are actu- 
ally cospatial. As shown in Fig. [T] a is the separation of two 
stellar centres, and D = 2.63 R0 is the distance between the 

two surfaces along the joint line of the stars. 

Because of the shor t orbit al period, IStruve fc Zebergsl 

(|l959l ) and lOian et all (|2002T l suggested that YY Gem 
should be in synchronized rotation. Therefore, we assume 
that the angular velocities of both stars have the same value 
as that of the co-rotating reference system. This means that 
the two stars rotate almost face to face except the differ- 
ential rotation. As a result of the differential rotation, per- 
sistent shear motion is imposed to the coronal loops whose 
footpoints are rooted at different latitudes. 

The MHD processes driven by the differential rotation 
are three-dimensional in nature. The addition of the two 
spherical boundaries makes the problem more unfeasible to 
solve numerically. To simplify the problem without losing 



the essence of physics, the spherical boundaries are treated 
as two parallel planes, which are tangent to the spheres at 
the equators as illustrated with the dotted lines in Fig. [T] 
The area of interest is ^ x ^ D and —2.5 Rq ^ y ^ 2.5 
Rq, where the s-axis is along the joint line of the binary, 
the y-axis is parallel to the axis of the rotation, and the 
z-axis is perpendicular to the x — y plane, with the origin 
of the coordinate being located on the equator of Mi star. 
Note that the vertical size is set to be larger than the stellar 
radius in order to minimize the influence of the top bound- 
ary on the numerical results. Moreover, all quantities are 
approximated to be invariant along the z-axis, so that the 
problem becomes 2.5-dimensional. In reality, the interface 
between the magnetic systems of the two stars has a limited 
extension in the z-direction. Therefore, it should be kept in 
mind that any feature in our numerical results should have 
a limited size in the z-direction. 



2.2 MHD equations 

We perform 2.5-dimensional (d/dz — 0) numerical sim- 
ulations in the Cartesian coordinates as descr ibed above. 
The M HD equations are slightly modified from lChen et al.l 
in order to include the gravity from both stars and 
the centrifugal force in the co-rotating system. For simplic- 
ity, heat conduction is not included. The resulting dimen- 
sionless MHD equations, which are shown belo w, are nu- 
merically solved by a m ultistep implicit scheme (|Hulll989l ; 
IChen. Fang fc HrtoOOT ). 
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where p, v x , v y , v z , tp, B z and T correspond to the dimen- 
sionless density, three components of velocity, magnetic flux 
function, z-component of the magnetic strength, and tem- 
perature, respectively, and po, t»o, ^0, -Bo and To are the 
characteristic values of the corresponding parameters, which 
are used to nondimensionalize the equations. The charac- 
teristic velocity vo is defined as the isothermal sound speed 
vq — \/RTo, where R is the gas constant for the fully ionised 
hydrogen and the length scale Lq is set to be 0.5 Rq. The 
dimensionless size of the simulation box is then ^ \y\ ^ 5 
and ^ x ^ 5.26. Since the flux function and the magnetic 
field are related by 



V x (tPe z ) + B z 



(8) 
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Figure 1. Sketch of the binary system. For the case of YY Gem, the two stars are identical. Here, Sx is the deviation of the x-coordinate 
from the centre of mass, n is the distance to the centre of Mi, and T2 is the distance to the centre of M2. The simulation box is indicated 
by the dashed box. 



where e z is the unit vector in the z-axis, we have tpo = BqLq. 

In equation g is the total acceleration of gravity 
contributed by both stars, 



9 = 



G M 

I _2 



Mi 



VnL, 







(9) 



where Mi and M2 are the masses of the two stars, n and 
T2 are indicated in Fig. [T] (note that such an expression is 
an extension of the gravity along the x-axis, which is valid 
since we are interested in the processes near the z-axis). In 
the same equation, F c represents the centrifugal force, 

G(Mi + M 2 )L 2 



F c = 



-5x, 



(10) 



where a is the distance between the two stars, 8x is the 
deviation of the ^-coordinate from the centre of mass. The 
plasma beta, /3o, and the magnetic Reynolds number, Rm, 
are expressed as 



0o 
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where fio is the magnetic p ermeability and T) is t he anoma- 
lous resistivity. Similarly to lChen fc Pr iest (2006), r\ is cho- 
sen to be a function of the current density, j, 



77omin(l, j/j c 
0, 
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(13) 



where j c is the critical current density. 

The time unit used in the description of the numerical 
results is the Alfven time-scale, ta = Lq/va, where va is 
Alfven speed, i.e., 



VA = 



(14) 




Figure 2. Initial configuration of the magnetic field (solid lines) 
and the distribution of the dimensionless density (gray scale). 



As iTsikoudi fc Kellettl (|200d ) pointed out, the corona 
of YY Gem is very hot, To can be taken to be sev- 
eral million kelvin. The number density (n = p/w p ) 
of the stellar corona is in the range of 10 14 - 10 18 
m~ 3 (iMonsignori-Fossi fc Landinil Il994l iMewe et al.1 [l995l; 
Schmitt et all 1 19961 ; ISchriiver fc Zwaanl I 2000l : iNess et al.1 
2002T I. The typical strength of the general magnetic field near 
the stellar surface is tens of Gauss. Therefore, in our stan- 
dard case, which is called case A2 later, the following charac- 
teristic values are adopted: no = 10 15 m~ 3 , To = 1.9 x 10 B K, 
/?o=0.01 (hence Bo=36 G), and j c = 10. The corresponding 
time-scale is ta = 139 s. 



2.3 Initial and boundary conditions 

Initially, the isothermal plasma is at rest, i.e., v x — v y — 
v z — 0, and T = 1. The density distribution is obtained by 
numerically calculating the hydrostatic equation 
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p ax 



(15) 



The initial magnetic field in this work is assumed to be 
potential, which is produced by two parallel line currents 
lying below the two surfaces, i.e., 



i> = - ln(V \x - xi) 2 + y 2 ) - ln(^(x-x 2 ) 2 +y 2 ), 



(16) 



where x\ and X2 are the ^-coordinates of these two currents. 
We take x\ = — 1 and x 2 ~ D/Lo + 1 in our work. The result- 
ing density distribution and the magnetic configuration are 
plotted in Fig. [2] where the magnetic f ield is similar to the 
configuration adopted in lFerre ira (1998). The magnetic con- 
figuration is taken in such a way that a magnetic null point 
exists between the antiparallel magnetic loops. So, when- 
ever the closed loops bulge, a current sheet forms natually, 
and magnetic reconnection is ready to commence (if the po- 
larity of one magnetic system expressed in equation (|16|l is 
flipped, the two loop systems become to repel each other, 
and no reconnection can happen). Similar magnetic connec- 
tivity with a nu ll point was also computed for the active 
RS C Vn system (lUchida fc Sakurai|[l983l ; lUchida fc Sakurail 
ll985l ; lBeaslev fc GiidellbOOOl ). 

Owing to the symmetry, calculations are done only in 
the upper half of the simulation area, i.e., ^ x ^ D and 
^ y ^ 2.5 R,0. The calculation domain is discretised into 
501 uniform grid points along the a>axis and 141 nonuniform 
grid points along the y-axis, with more grid points concen- 
trated near the a-axis. The top side of the simulation box 
(y — 5) is treated as an open boundary, and the bottom one 
(y = 0) is a symmetry boundary. The right and left, which 
correspond to the surface of the binary stars, are line-tied 
boundaries, where all quantities are fixed except that a typ- 
ical differential rotation is imposed at each surface in the 
following form: 



I left 



'right 



SIR, 



sin 2 i + SIR, i < ir/2, 
i > tt/2, 



OR(l-asin 2 i)\JT 
—SIR, 



SIR, i < tt/2, 
i > tt/2, 




where SI is the angular velocity of the binary, a is the dif- 
ferential rotation rate, and i, the latitude, is related to the 
coordinate y by y = Rsini. Beyond y = R, v z is set to be 
the value at y — R. In this paper, a is set to be 0.5. As we 
will discuss later, the value of a does not affect the main 
results in our numerical simulations. 



3 NUMERICAL RESULTS 
3.1 General evolution 

The distributions of the magnetic field, the temperature, 
and the velocity field at several selected times in case A2 
are shown in Fig. [3] As the shear motion is imposed at the 
left and right boundaries, footpoints of the coronal loops 
are dragged to move in the z-direction, which increases the 
local magnetic pressure. As a result, both the left and the 
right magnetic loop systems begin to inflate and approach 
each other due to the increasing pressures near the roots of 
the loop systems (see the snapshot at t — 83 ta)- Such an 
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Figure 4. Temporal evolution of the reconnection rate in case 
A2, where R r is the reconnection rate, and ta = 139 s. 



inflation of magneti c loops due to shear motio ns was demon- 
strated on the Sun (|Barnes fc Sturrock|[l972l ). The two loop 
systems collide near the magnetic null point, by which a 
current sheet forms. When the increasing current density 
exceeds the prescribed threshold (j c ), anomalous resistivity 
is excited, which triggers the occurrence of magnetic recon- 
nection between the two loop systems (see the snapshot at 
t = 109 ta). The temporal evolution of the reconnection 
rate, R r , which is defined as dip/dt at the magnetic null 
point, is shown in Fig. [4] It is seen that the reconnection 
rate presents a periodic behaviour. With the wavelet anal- 
ysis, the period is found to be 54 ± 6 min. Note that, with 
the 2.5D approximation in our model, the reconnection be- 
tween the two closed magnetic systems, which are rooted on 
the surfaces of the binary stars, will go on as long as the 
differential rotation continues. In reality, the differential ro- 
tation of the two spherical stellar surfaces would take the 
two closed magnetic systems apart, which would cease the 

(17) reconnection after a certain time. The 2.5D model is valid 
here since we are interested in the magnetic reconnection 
only during first several periods. 

(18) In order to compare th e numerical result with the U- 
band observations made bv lDovle et alj (|l990h . we synthe- 
sised the UV emission in the wave- band of 3650 ± 10A 
from our numeri cal results with the help of the CHI- 
ANTI database jDere et al.l Il997l ; lLandi et al] 120061 ; see 
IChen fc Priest! 120061 for details'). The calculated UV right- 
curve of the whole simulation area is plotted in the left panel 
of Fig. 5, and its wavelet spectrum is presented in the right 
panel. Note that since the reconnection rate and the UV 
emission both are small in the first 60 ta, we displayed the 
result after t = 60ta- It is seen that the UV flux does show 
a periodicity, and the period is centred at 52 min, very close 
to that of the reconnection rate, i.e., 54 min. 



3.2 Parameter survey 

In order to investigate how the choice of the characteris- 
tic values of various parameters may affect the results, we 
have made numerical simulations with different characteris- 
tic values, which are listed in Table [1] In cases A1-A5, only 
To is changed. Similarly, cases B1-B5 are for /3o dependence, 
cases C1-C3 for po, and cases D1-D3 for j c , respectively. In 
all cases, Lq is fixed to be 0.5 Rq. Numerical simulations 
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Figure 3. Magnetic field lines (solid lines), velocities (arrows), and temperatures (colour) at five selected times in case A2. 
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Figure 5. Left panel: Temporal evolution of the flux in the [/-band , which is derived by using the CHIANTI database. Right panel: 
The wavelet spectrum corresponding to the [/-band flux, which has a period centred around about 52 min. 
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Table 1. Parameters used in different simulation cases. 



Models 


T (xlO 6 K) 


A) 


lgn /m 3 


jc 


P (min) 


Al 


1.5 


0.01 


15 


10.0 


70 ± 9 


A2 


1.9 


0.01 


15 


10.0 


54 ± 6 


A3 


2.5 


0.01 


15 


10.0 


41 ± 4 


4 1 


6 . U 


n m 


la 


inn 


n rr | o 


AO 


l.U 


n m 


"1 p, 
10 


inn 


ol it O 


Bl 


2.0 


0.05 


15 


10.0 


96 ± 15 


132 


2.0 


0.025 


15 


10.0 


78 ± 10 


B3 


2.0 


0.01 


15 


10.0 


52 ± 6 


B4 


2.0 


0.0075 


15 


10.0 


47 ± 4 


B5 


2.0 


0.005 


15 


10.0 


38 ± 4 


CI 


2.0 


0.01 


14 


10.0 


52 ± 6 


C2 


2.0 


0.01 


15 


10.0 


52 ± 6 


C3 


2.0 


0.01 


16 


10.0 


52 ± 6 


Dl 


2.0 


0.01 


15 


1.0 


52 ± 6 


D2 


2.0 


0.01 


15 


5.0 


52 ± 6 


D3 


2.0 


0.01 


15 


10.0 


52 ± 6 


indicate that, similar to 


case A2 


as described 


in the 


previ- 



ous subsection, all other cases present periodic behaviours. 
The corresponding period of each case is displayed in the 
last column of Table [T] It is found that, with To, /3o, po, 
and j c as independent parameters, the period of the flaring 
is strongly dependent on To and /3o, and an increasing T 
or decreasing /3q results in a decreasing period. When po in- 
creases, the magnetic field Bo increases as well in order to 
keep /3o constant (correspondingly, the Alfven speed is con- 
stant). As a result, the period of the reconnection rate does 
not change. It is also seen from Table[T]that j c has little effect 
on the period, which implies that the magnetic reconnection 
is modulated by some external process. 

Besides the four parameters mentioned above, we also 
examine other parameters related to the binary itself. In 
case A2, by changing the values of a (a = 0.2, 0.3, 0.4, 

0. 5, 0.6, 0.7, 0.8), we find that the results keep invariant, 

1. e., P = 54 ± 6 min. In other words, the flaring period has 
nothing to do with the speed of the differential rotation. 
As mention ed above, the mass and the siz e of YY Gem are 
taken from iBrancewicz fc Dworakl lll980j). However, other 
authors derived slightly different values ( Leung fc Schneider! 
ll97Sl;IChabrier fc Baraffelll995l : iTorres fc Ribaj|2002l ). e.g.. 



in iLeung fc Schneider! (| 19781 ). Mi = 0.62 M , M 2 = 0.57 
Mq, Ri = 0.66 R,0, i?2 = 0.58 R©, the separation of two 
components a — 3.9 Rq, and the orbital period P = 0.82 
d. Through test calculations, we find little change in the 
results. 



4 DISCUSSIONS 

4.1 What modulates the magnetic reconnection? 

In order to understand how the magnetic reconnection is 
modulated, the temporal evolution of v x along the z-axis 
is shown in Fig. [6] where the value of v x is indicated by 
the grey scales. As the reconnection rate R r is equivalent 
to v x B y along the x-axis, it is expected to see that v x at 
a given point also shows an oscillatory evolution, with the 
same period as R r , i.e., 54 min. Furthermore, a wave pattern 
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Figure 6. Temporal evolution of v x along the x-axis. 
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Figure 7. Temporal evolution of reconnection rate in the case 
that the horizontal size of the simulation box is cut down. 



is discerned in Fig. U5] where the waves are seen to travel back 
and forth in the a;-direction, bounded by the surfaces of the 
two stars. Since the waves propagate perpendicular to the 
magnetic field lines, we postulate that they are fast-mode 
magnetoacoustic waves, and that the periodic magnetic re- 
connection in the numerical results is modulated by such 
waves. 

To verify that the magnetic reconnection is modulated 
by magnetoacoustic waves bouncing back and forth between 
the two stellar surfaces, we shift the right boundary of case 
A2 leftward and the left boundary rightward by 0.8 simul- 
taneously, by which the separation between the two stellar 
surfaces is cut down from 5.26 to 3.66, i.e., a decrease with 
a ratio of 30 per cent. The corresponding evolution of the 
reconnection rate is plotted in Fig. [7] Note that the posi- 
tive R r means that the field lines in the two loop systems 
are reconnecting, whereas the negative R r means that the 
transverse field lines are reconnecting. 

It is found that the resulting period of magnetic recon- 
nection becomes 34 ± 6 min. This means that the period 
decreases by 37 per cent, a ratio close to that of the de- 
crease of the separation between the two stellar surfaces. 
This implies that the magnetic reconnection is modulated 
by waves trapped between the two stellar surfaces. 

To verify that the travelling waves are fast-mode waves, 
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Figure 8. Temporal evolution of the wave front. Only the right 
half is plotted owing to the symmetry of the YY Gem system. 

here we calculate the trajectory of fast-mode waves propa- 
gating from the reconnection point (x = 2.63), through the 
right boundary (x = 5.26) and then back to the reconnection 
point along the x-axis. For a given number n, the travelling 
time (T) and the position of the wave front (a;) are deter- 
mined by the following equations 

n n 

T = y^Ax/vf, x- 2.63 = Ax, (19) 

where Ax is the uniform grid spacing in the 2>direction, and 
Vf is the fast-mode wave velocity. The procedure is described 
as follows: Starting from the reconnection point (x=2.63), 
the wave front propagates rightward by Ax each step, the 
time interval At is equal to Ax divided by Vf, so the total 
travelling time is the sum of At, while the position is the 
sum of Ax. Note that since both the magnetic field and the 
plasma temperature are changing with time, we update the 
magnetic and thermal parameters every 1 ta when calculat- 
ing Vf. Besides, when the wave reaches the right boundary 
and then bounces back, Aa; takes a negative value. After 
the wave front returns to the reconnection point, the wave 
should propagate to the left half area of the simulation box. 
Considering the symmetry of the problem, i.e., there is a 
symmetric wave propagating rightward from the left half 
area, we repeat the above calculation in the right half re- 
gion. 

Figure. [5] shows the temporal evolution of the position 
of the wave front in the right half region only. It is found 
that the trajectory of the fast-mode wave, from the recon- 
nection point to the right (or left) boundary and then back 
to the reconnection point, is quasi-periodic, with a period 
of ~ 23 ± 2 ta, i.e., ~ 53 ± 5 min, very close to the period 
of the magnetic reconnection rate presented in Section 3. 
Therefore, we are convinced that the magnetic reconnection 
is modulated by fast-mode waves. The underlying process 
can be understood as follows: When symmetric waves, intro- 
duced either by the shear motions of the stellar surfaces or 
by the magnetic reconnection, collide near the reconnection 
point, the reconnection rate is enhanced to form a peak. Af- 
ter the collision, the two waves propagate to the right and 
left boundaries and are bounced back, respectively. When 
they collide again near the reconnection point, the recon- 
nection rate is enhanced again. Therefore, the period of the 



magnetic reconnection is determined by the travelling time 
of fast-mode waves from the reconnection point to the right 
(or left) boundary and then back to the reconnection point. 

4.2 Empirical formula for the period 

As shown in Section 3, among the selected four parameters, 
the period of the reconnection rate (P) is explicitly depen- 
dent on To and /3o- The underlying process can be under- 
stood as follows: When To increases while keeping /3q and 
po constant, both the sound speed and the Alfven speed in- 
crease, which leads to a larger fast-mode wave speed. As a 
result, the period of the reconnection rate decreases; When 
Po decreases while keeping To and po constant, the sound 
speed does not change, while the Alfven speed increases, 
which also leads to a larger fast-mode wave speed and a 
shorter period; When po increases while keeping Po and To 
constant, both the sound and the Alfven speed do not change 
(the Alfven speed does not change since the magnetic field 
Bo increases). As a result, the period keeps constant. 

In order to derive an empirical formula to relate P to 
To and /3a, we simulate a series of cases with different To 
and Po. The corresponding variations of P with To and /3o 
are displayed by the squares in Fig. [9] It is seen that lg P is 
almost a linear function of lg To and lg Po , which means that 
P scales with To and Po as P oc Tq 71 /?^ 2 . By fitting the data 
points with the solid lines in Fig. [9] we obtain 71 = —0.91 
and 72 = 0.42. Therefore, we have P oc T "°- 91 /3g 42 , i.e., 

P~ pl A2 T ' 0A9 B^ &i . (20) 

According to our postulation that the magnetic re- 
connection is modulated by fast-mode waves, the period 
P should be proportional to 1/vf, where = <? s + v\, 
c? s — 7-RT0, and v\ = Bof(fiopo)- In the extreme case of 
zero magnetic field, we have P oc T -0 ' 5 ; In the case of ex- 
tremely strong magnetic field, we have P oc pJP-Bq" 1 - The 
slight deviation of the power indices in our numerical results 
from those in the two extreme cases further implies that it is 
the fast-mode waves that modulate the magnetic reconnec- 
tion. Based on this scaling law, the observed period of the 
flaring rate in YY Gem, i.e., 48 min, corresponds to a char- 
acteristic magnetic field of 41 G, which is the value of the 
magnetic field in low corona, providing the number density 
of the plasma in the inter-binary space is about no ~ 10 15 
m" 3 . 

4.3 Other points 

The period ic flari ng of YY Gem, which w as found by 
iDovle et al.l (|l990h . was not confirmed in IStelzer et al.l 
i|2002l ). The discrepancy might be attributed to two factors: 
(1) Probably the two observations were made at different 
phases of the stellar magnetic cycle. In order to have the 
large-scale magnetic reconnection in the inter-binary space, 
the magnetic configuration of each star should be close to 
the bipolar type, with interconnecting magnetic field lines 
linking the polar regions of the binary stars, as shown by Fig. 
[2] Such a magnetic configuration often appears near activity 
minim um; (2) The periodic flares observed by IDovle et al.l 
( 1990) and confirmed by our numerical simulations, are pro- 
duced in the inter-binary space, therefore, are large-scale in 
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Figure 9. Upper panel: Relation between the period (P) and the 
temperature (T); Lower panel: Relation between the period (P) 
and /3q . The solid lines are the corresponding linear fitting of the 
data points. 



nature. Without doubt, there should be small-scale flares 
that occur near the stellar surface and are due to the re- 
connection of the magnetic fields of a n individual star. The 
flares observed bv lStelzer et al.l (|2002h probably correspond 
to these small-scale flares. 



5 CONCLUSIONS 

Observations indicated that there may exist periodicity in 
the occurrence of strong flares from YY Gem. Using 2.5- 
dimensional MHD numerical simulations, we confirmed the 
periodicity of the flaring in the inter-binary space of YY 
Gem. On the basis of the analysis of the simulation results, 
the following conclusions can be drawn: 

1. Owing to the differential rotation of the stellar sur- 
face, coronal loops in each star inflate and approach each 
other. The resulting magnetic reconnection generates the 
inter-binary flares, which are large-scale in nature. Fast mag- 
netoacoustic waves, which are trapped in the inter-binary 
space between the surfaces of the two binary components, 
modulate the magnetic reconnection, producing the periodic 
behaviour of the flaring rate. 

2. The period scales with the typical coronal density 
(po), temperature (To), and the coronal magnetic field (Bo) 

as p~ Po ,A2 t- OA9 b-°- s \ 
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